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Abstract. In this paper we applied the contour integral method for the zeta function 
associated with a differential operator to the Laplacian on a surface of revolution. Using 
the WKB expansion, we calculated the residues and values of the zeta function at several 
important points. The results agree with those obtained from the heat kernel expansion. We 
also obtained a closed form formula for the determinant of the Laplacian on such a surface. 
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1. Introduction 



Spectral functions, namely functions of the spectrum of, usually, second order differential 
operators play an important role in mathematics and physics. The possibly most important 
such spectral function is the zeta function which directly relates to topics such as analytic 
torsion [IQ], the heat kernel [271 S3], Casimir energy (U El H31 EH [36], effective actions and 
quantum tunneling jZJ |TTl El H51 ESI I2H 123 S2] and the study of critical metrics [TOj 1381 139] . 
The relationship between the zeta function and all these topics is established by the fact that 
very specific properties of the zeta function encode the information needed in these areas. 
For example residues and values at specific points contain information about the short time 
rn ! heat kernel asymptotics with all the geometric information it contains, the (finite part of 

its) value at s = —1/2 relates to the Casimir energy, and its derivative at s = is relevant 
in all the other subjects mentioned above; see, e.g., |29j. It is therefore no surprise that an 
enormous effort is made to understand spectral zeta functions as completely as possible. 

The situation is particularly satisfying for the one- dimensional case where functional de- 
terminants can be obtained in closed form as boundary values of a suitable initial value 
problem El [T71 [22], [231 [261 US]. A contour integral approach established in [311 [32] was 
used to rederive this case and generalize it to arbitrary boundary conditions. This approach 
has the advantage that it allows for the systematic study of higher dimensional configura- 
tions as long as the Laplacian separates in a suitable fashion. Examples are Laplacians with 



spherically symmetric potentials [TS1 [T9] and Laplacians for spherical suspensions [24|. In 
those cases the analysis of the associated zeta function is based on explicitly known special 
functions, namely Bessel functions and Legendre functions. 

Another example, studied in the present paper, is the Laplacian on a compact surface of 
revolution. These surfaces are geometrically interesting enough to exhibit many features that 
appear in the study of spectral functions and their asymptotics, while at the same time the 
rotational symmetry allows the explicit calculations that are necessary in order to understand 
meromorphic extensions of the zeta function to regions outside of the domain of convergence. 
The complication with this example is that no known special function emerges as a solution 
to the relevant initial value problem because the function / whose graph is used to generate 
the surface of revolution is kept fairly general. As a result needed asymptotic behavior has 
to be determined from scratch starting from an ordinary differential equation. The WKB 

method is our method of choice [31 [35] and it allows us to obtain the values or residues of the 
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zeta function at several important points, in addition to the determinant of the Laplacian. 
The relation between the values of the zeta function and the global geometrical properties 
of the surface is revealed. Furthermore it is notable that the result for the determinant is a 
closed-form formula. 

The structure of the paper is as follows. In Section [2] we present the geometric prelimi- 
naries, define the generalized zeta function, and reformulate it as sum of contour integrals. 
The zeta function is written as the sum of a one-dimensional term and a two-dimensional 
term, which are calculated in Section [3] and Section H] respectively. The values and residues 
of the zeta function at several selected points are calculated along with the determinant of 
the Laplacian in Section |H We conclude in Section with a comparison of our method and 
other approaches. 



2. Zeta function for a surface of revolution 



2.1. Laplacian on a surface of revolution. We consider the Laplacian on a surface of 
revolution with Dirichlet boundary values. Let / G C 2 (a, b) be a positive function from [a, b] 
to M, and let M be the compact surface with boundary that is generated by revolving the 
graph of / around the x-axis. With the metric induced by the Euclidean metric of M 3 , (M, g) 
becomes a Riemannian manifold. In the natural coordinates (x, 9), the metric is 



1 + /'(a;) 2 








Recalling the formula for the Christoffel symbols in terms of the metric, 



1 2 

9 5^ ' 



M 



dgu dg^ 
dxi dxi 



dgij 
dxi 



we have 



r 1 



f'f 



r 2 

1 12 



r 2 

1 21 



ii 1 + fa 

and the others are zero. We then calculate that 
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The eigenvalue equation is Am = —Am. The following facts are standard [25j 

(1) The eigenvalues are real and positive; if labeled and ordered as < Ai < A 2 < . . ., 
then Afc — > oo. 

(2) For each A&, the corresponding eigenspace in L 2 (M) is finite-dimensional. 

(3) The spectrum of M is discrete. 

We remark here that zero is not an eigenvalue of the Laplacian if Dirichlet boundary values 
are taken. 



2.2. Generalized zeta function. Our spectral function for the Laplacian generalizes the 
well-known Riemann zeta function, 



oo 

Us) = 

k=i 



(2) 
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which converges for dt(s) > 1. For the Laplacian on a surface of revolution with Dirichlet 
boundary condition, if the eigenvalues are labeled and ordered as < A x < A 2 < • • ., the 
generalized zeta function is defined to be 

c(*)= E v, ( 3 ) 

A fe GSpec(A) 

for those values of s £ C for which this sum converges. Weyl's asymptotic formula can be 
used to determine the domain of convergence of Eq. ([3]) [12]. On a compact Riemannian 
manifold M of dimension n, 

2 2"7r5r(^ + l) 

From this we determine that the zeta function is defined for complex values of s having 
3fJ(s) > n/2, in our case 9ft(s) > 1. However, just as for the Riemann zeta function, much 
information of geometric and physical significance is contained in £(s) for complex values of s 
outside of this region, ((s) has a simple pole at s — 1. The holomorphic continuation of ((s) 
to 3?(s) < 1 turns out to have simple poles at s = | — n, n £ N. Here we obtain expressions 
for the residues of ((s) at s = 1, |, — |, and its value at s — 0, in terms of geometrical 
properties of the surface. Furthermore we compute the determinant of the Laplacian using 
the relation with the zeta function given by [161 EE, HQ] 

lndet A = -C'(0). (4) 

We obtain a formula for C'(0) given entirely in terms of the function / whose graph generates 
the surface of revolution. 



2.3. Reformulation of the zeta function as a contour integral. A fundamental diffi- 
culty in studying generalized zeta functions is that the eigenvalues are not known explicitly 
except for special cases. However, as described in Ref. [2S1 ED], the zeta function can be 
reformulated as a contour integral, which is highly suitable to such investigations. Here we 
adapt this method to the Laplacian on a surface of revolution. 

The rotational symmetry of M suggests the use of separation of variables to describe the 
eigenvalues of the Laplacian. The eigenfunctions of the Laplacian on the cross-section S 1 , 
denoted by Uk(9), k £ Z, are e tk0 , with eigenvalue —k 2 . Let u(x, 6) = (p{x)uk{0). By Eq. (JTJ, 
u is an eigenfunction of A on M with eigenvalue A if and only if is a nontrivial solution to 

fix) + ( j " + (a - |T) (! + f' 2 )^) = °> (5) 

satisfying the boundary conditions 

0(a) = 0, 0(6) = 0. 

Denote by 0fc(A; x) the solution to the initial value problem consisting of the same differential 
equation (jSJ), but with initial conditions 

0(a) = 0, 0'(a) = 1. (6) 

Then for each integer k, A £ C is an eigenvalue if and only if 0fc(A; b) = 0. For each fixed 
integer k, order the eigenvalues and denote them by < A^i < A^ < — Since 0fc(A; b), 
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regarded as a function of A, vanishes at and only at \k, n , n £ its canonical product 
representation [T] is 

00 / A 
^(A;6)=0 fc (O;6)TT 1- — 

n=l V Xk ' 

The convergence of the infinite product is guaranteed by Weyl's asymptotic formula. Let 
By the Residue Theorem, we may express ((s) as the sum of contour integrals. 

J1 k=l Jl 

where Y is any counterclockwise contour that encloses all the Afc jn 's but not the origin. The 
factor 2 in the second term above comes from the fact that Afc, n = \-k, n - Following Ref. [30] 
(Fig. 3), the contour can be chosen to go around the negative axis, from —00 to above 
the real axis, along which A = e l7T y, and from to —00 below the real axis, along which 
A = e~ tn y. Using this contour, Eq. (J7J) is simplified to be 

C( s ) = Ci(s) + C 2 (s), (8) 

where 

shiTrg f°° _ s d\n D (-y) 

CiOO = / y — j dy, ( 9 ) 

Jo dy 

Note, that Eq. (|T0|) should be understood in the sense that the analytical continuation of 
the integral is constructed first and the summation of that result is performed afterwards. 

Next, in Section[3], we will use the WKB method to solve Eq. (jSJ) with k = 0, and compute 
important values of Ci( s ) an d Ci(0)- I n Section H] we will adapt the method to the series over 
k, and compute (2(3) along with £2(0)- 

3. Zeta function and derivative with k = 

By Weyl's asymptotic formula Eq. Q converges for 3?s > ~. To obtain the meromorphic 
extension of (j(s) to 9Rs < |, we apply the WKB method to determine the asymptotics of 
0o(— y; b); see for example [21 El] for a discussion of the WKB method. First, we make the 
change of variables y = z 2 in the differential equation (jSJ) with k = 0. It becomes 



*) + - T~T~co - + Z' 2 )^(^) = 0- (11) 



//, / £ J J \ . // .. i 2, 1 . 

/ 1 + f 2 , 

Denote the coefficient of </>'(x) in Eq. ( |TT1) by u(x), 

r ff" 



U\X) 



f 1 + f 2 ' 

Let 



<f>(x) = exp ( / a) - ^}da ) . (12) 
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Substituted into Eq. (TTTT) . the equation for S(z;x) is 

S>(z;x) + S(z;xf = (l + nz' + ^ + ^, 

where the prime indicates differentiation with respect to x. The solution can be written as 
asymptotic expansions in z. It has two branches, namely, 



S±(z; x) = ±s_ lZ + s ± — + . . . = ±(s_iz + — + ...) + (s + — + .. (13) 

z z Z A 

where the coefficients s n (x) can be computed recursively. The first three coefficients are 

f'f" 1 f' 2 , 1 /" 



/2 



SO 



Si 



2(1 + /'*)' ^ 4/(l + /T' 

Applying the initial condition 0(a) = and </>'( a ) = 1, the solution to Eq. (ITT]) is 



(14) 



i r* , I, e ^ s+{ - z '' a ' )do ' - e^ s -^ da 



S+(z; a) — S-(z; a) 

I — e ~ fa[S+{z;°-)-S-(z;a)]dcT 



S + (z; a) — S^{Z] a) 



S2 



Let S± = ±Si + S 2 , where 

S^z; x) = + — H , 5 2 (z; x) = s + ^ + ■ ■ 

Notice that 

S; + ^ = ^ + 5 2 ^^ + 25 1 5 2 = 0^1n5 1 (^;a)-ln5 1 (z;6) =2 / S 2 (z;x)dx. 

J a 

(f) (—z 2 ;b) can be written as 



(15) 



cj) {-z 2 ;b) = e/. 6 [Si(*:*)-s«(*)]^ 



2 A /S , i(^; a) Si (2; 6) 



(16) 



Now we apply Eq. f[T6l) to the contour integral Eq. ([2]). Noticing that Eq. (lT3"j) is invalid as 
z — >■ 0, we write the zeta function as 



Ci(s) 



sin tcs 



71 



1 _ s d\nP (-y) ^ 

y ; dy + 



dy 



00 _ 2s d\nD (-z 2 ) 



dz 



By Eq. (JT6 



InZU-z 2 



6 r 



Si(z; x) - 



u\x) 



dx 



dz 



In Si (z; a) + In S\ (z; b) 



(17) 



+ ln 



X _ e -2 f*Si(z;x)dx 



ln(20 o (O;6)). 



To compute dD (—z 2 )/dz, we need the asymptotic expansion of dS\(z; x)/dz, which can be 
derived from Eq. (fl5l) . and the expansion of dhiSi(z;x)/dz. Let 



<91n Si (z; x) 9 
<9z <9z 



ln(s_iz) + 



Sl 



S_lZ z 



+ 



t (x) h(x) 
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where 

t (x) = 1, tx(x) 
Substituted into Eq. (TTST) . 



-2 



si 



./ 



/2 



r 



S_! 4/2(l + /«) 2/(l + /«)2' 



(19) 



lnDnf-z 2 



s_i(x)<ix 



Si(x)da; 



l t (o) +t (6) _ l ti(o) + ti(6) 
2 



2^ ^ z J 

Now we substitute the equation above into the zeta function in Eq. (j!7p . The last integral 
above is uniformly bounded for s with 9fts bounded below, while the integration from to 
1 in Eq. (ITTj) is uniformly bounded for s bounded above by 1 — e for any e > 0. Using the 
meromorphic extension of the following integrals to s G C, 

- 2s z- n dz=- r , n = 0,1,2,... (20) 



2s + n-r 111 

we can compute the residue of the integral from 1 to 00 in Eq. (JIT]) at s = |, 0, — |, etc. If 
s in an half integer, sin7rs ^ 0, and Ci has a single pole at s for a general /(a;). If s is an 
integer, shi7rs has a single zero at s, and so £i has a finite value at s. For example, using 
Eq. ( 120]) with n = 0, 1, 2, 3 respectively, and plugging into Eqs. ({14]) and (jT9]) . we have 
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2f(fc) 





f'(x) 2 dx, 



f" 



4/(1 + 7/2)3/2 

2/"(a) 
/(a)(l + /'(a) 2 ) 2 



fix, 



(21) 
(22) 
(23) 



+ /(6) 2 (1 + /'(&)') /(&)(1 + /'W) 2 . ' (24) 
The calculation can be continued for s < — 1, but the expressions become very complicated 
and we do not include them here. 

Derivative of Ci( s ) at s = 

To compute ([(0), we use integration by parts. Rewrite Ci( s ) as 



&0O 



sin tys 

s I y 
~ Jo 

sm7rs '"' 



7T 



-a-1 



In D (-y)dy 
s - 1 \nD (-y)dy + 2 



w -2b-1 



(25) 



y " * mu {—y)ay + 1 1 z ~ - \nD (—z )dz . 
'0 Ji 

Since \nDo(—y) = 0(y) for y — > 0, the integration from to 1 is bounded for s = 0. For the 
integration from 1 to 00, by the asymptotic expansion Eq. ( |T5l) . 

Si(z; x) = s-i{x)z + Oiz- 1 ) = s- 1 (x)z ■ [1 + 0{z~ 2 )l 
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as z — > oo. Substituted into Eq. ([IB]) . 



\nD (-z 2 ) = 
-0(z- 2 )+\n 



" r s _ 1 ( x )^ + 0(z- 1 )-^ 
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Substituted into Eq. fl25l) . as s — >■ 0, 



2s- 



2s 2 
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2l 
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ln(20 o (O;6)) 



This gives Ci(0) = — |' as * n (^P . and furthermore the following theorem. 
Theorem 3.1. The derivative of the zeta function in the one- dimensional case is given by 

ln/(a)+ln/(6) 
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Proof. By the calculations above, we have 



Cl(o) 



u(x)dx + 



In s_i(a) + In s_i(fo) 



+ ln(20 o (O;6)) . 



io(0; b) can be computed exactly, 



0o(O;6) 



'1 + /' (a^ 2 



where A is defined by Eq. (|2Tj) . Combined with u(x) = (In /-In v/1 + J*)', Eq. (j2BJ) follows 
from a simplification. □ 
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4. Full zeta function and derivative 

For k 7^ 0, we make the change of variables y = k 2 z. The differential equation for (j){x) 
becomes 



<f>"{x) + u{x)4>\x) - k\l + f) z + 



P 



(x) = 0, 



(2? 



where u(x) is as defined in Sectional The zeta function in Eq. (TTUj) becomes 



Ca(«) 
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7T 
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00 _ 5 dlnD fc (-fc 2 ^) 

^ dz 



(29) 



Because the integrand depends on fc, the sum and the integral do not factor, but the sum 
nevertheless resembles the Riemann zeta function. This is a principal motivation for the 
change of variable. We define S(k, z; x) as in Eq. (TT2|) . and Eq. ffToT) holds with O (— z 2 ; b) 
replaced by (pk(—k 2 z; b). The equation for S(k, z\ x) is 



S'(k, z; x) + S(k, z; xY 



1 + / 



12 



and the solutions are written as asymptotic expansions in k, 



Si(k, z; x) = s-i(z; x)k H — — p h ■ 



(30) 



The first two coefficients are 



s-i{z;x 



where 



./ 
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It- 4 f 
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1 /" 
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Recall that D k (-k 2 z) 
In D k (-k 2 z) 

Similarly, let 

din Si(k, —k 2 z; x) 



t(x,z) = zf(x) 2 . 
^ k (-k 2 z;b)/M^b). By Eq. (EED, 

In S\(k, z; a) + In S\(k, z; b) 



(31) 
(32) 
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+ ln 



I _ fa Si(k,z;x)dx 
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;i + t) 4 i + f 2 4(i + t) 3 (i + / 
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(34) 
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Substituted into Eq. (I33l) . 
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dt 



dx + k~ l / f{x) 
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Substituted into Eq. (J2H 
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S\(k, z; x)dx. 



As before, Cr is the Riemann zeta function. Since Si (A;, z; x) = 0{k\fz) for large or z, the 
last term above does not contribute to the residues or values of (^{s) for s being integers or 
half integers below (including) 1. Using the following standard formulas [33] , 

1 



sin its 



7T = r( s )r(i- s ; 



and 



oo z - s 



we obtain 

Us) = U2s-1) 



r( s -| 
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5s + 1 f h fixf^f'ixf 
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-( R (2s + 2)f(a) 2sS - 
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4 Ll + /' 2 H 4 + (l + f 2 (a)) 
f 2 (6) 5s + 3 + fib) fib) 



-C fl (2 S + 2)/(6) 



2s ' 



_l + f 2 (6) 4 T (l + /'2 (6)) 2 

For s = 1, ~, 0, — |, etc, the residue or value of C 2 (s) may have two contributions from the 
series above, since both (r(s) and T{s) have poles. For example, using properties of the 
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Riemann zeta function, Cr(0) = ~\ an d Res Cr(1) = 1> we have 

ResC 2 (l) = ^/ f(x)^l + f'(x) 2 dx, 

f(a) + fib) 



ResC 2 (^) = ~ f y/1 + fWdx 

* ^ J a 



Ca(0) 



(35) 

(36) 
(37) 
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1 1 f'\a) 



f 
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1 /" 

+ T- 



pjr+r 4/(1 +/«)»/» 



+ 



/»/(a) 



32/(a) L 81 + f 2 (a) (l + /' 2 (a)) 2J 32/(6) L8 1 + /«(&) (1 + /*(&)) 



1 /"(&) /"(&)/(&) 



+ 



.(38) 



Adding up Ci( s ) an d C2(s), and recalling that Ci(l) is finite, we obtain the following theorem 
for the full zeta function. 



Theorem 4.1. For the full zeta function of the Laplacian on a surface of revolution, 

r-b 



ResC(l) = \j f(xWl + f>(x¥dx. 

/(a) + fib) 



ResC 



(39) 



(40) 



Res C {-\ 



32/(a) 



C(0) = 0. 

/"(a)/(a) 
l + /«(a) ' (1 + /'W 



1 f 2 (a 



+ 



32/(6) 



1 A&) 



+ 



(41) 
f"{b)fib) 



1 + /*(&) (l + f 2 (6)) 2 
(42) 



The geometric interpretation of the formulas above will be discussed in Section [51 
Derivative of £2(5) at s = 

The derivative of £2(3) at can be calculated in a similar but more subtle way. The 
following theorem is the main result of the paper. 

Theorem 4.2. The derivative at zero of the holomorphic extension of the zeta function on 
the two-dimensional surface of revolution is given by 



£W dx + i 

/(x)vTTTW 2 



f"(x) 



(l + /'(x) 2 ) 3 /2 

where A is defined in Eq. ( J271) . and <fi is the Euler function 0, [37] defined by 

00 

m = li(i-?*). 

k=l 



dx, (43) 



(44) 
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Proof. Integrating Eq. ( 129|) by parts, we find that 
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C2 ( s ) = 2s S -^J2 k ~ 2S L z- s - l \nD k {-k 2 z)dz 



k=l 



(45) 



Si(k, 0; x) and 0^(0; b) can be computed exactly, 



/(«) 



fix) xA+TW 2fc 

where A is defined by Eq. fl27|) . This allows us to write 
lnDtf-fc 2 *) = 



6 1 S 1 (k,z;a)S 1 (k,z;b) , . 1 - e -2/a*(M;*)* 

oi(k, -2; x)ax — — — In - - - — ; „ h - : ; + In ■ 



2 S , i(ib J 0;a)S , i(A:,0;6) 



By the asymptotic expansion in Eq. ( 1301) . we write 5i(/c, 2; x) as 

Si (z\ x) 

Si(/c, z; x) = s-i(z; x)k H ^ h 2; 2) 



1-e 



(46) 



where the remainder R(k,z;x) = 0(y~ 3 / 2 ) = 0(k~ 3 z~ 3 / 2 ) as y — > 00. In addition, since 
Si(k,0;x) = s_i(0;x)/c, we have si(z;x) = O(z) and R(k, z;x) = 0{k~ 3 z) as z — > 0. It 
implies that ln5i(fc, 2; x) can be written as 

In Si(k, z\ x) = ln(s_i(;z; x)k) + D(k, z\ x), 

where D(k,z;x) = 0(y _1 ) = 0{k~ 2 z~ l ) as y — > 00, and D(k,z;x) = 0{k~ 2 z) as z — > 0. 
Substituting Eq. f|475|) into Eq. f|45|) and using the expressions for S\(k, z\ x), we have 



„ , . sm7rs , 
( 2 (s) = 2s J2 k 

71 * — ' 



-2s 



k=l 



POO 








"1 


f / & " 


/ z ~ 




/Q 




[J a 



y/l + f 
f 



k(VtTT - 1) + + i2(Jfe, z; x) 



dx 



\n(zf(a) 2 + 1) +ln(^/(6) 2 + 1) D(k, z; a) + D(k, z;b) 1 - e~ 2 J*Si(M;*)«fc 

I In t tt ^ cLz . 



4 2 !_ e -2Afc 

The asymptotics of i2(/c, z; x) and z; x) as y — > 00 and as z — > implies that 



fc=i 



is analytic for — ^ < Sfts < 1. Therefore, 



2s- 



sm7rs 



7T 



2* 



fc=l 









r / b 








[7a 



/ 



-fc(V*TT- 1) + 



si(z;xj 

k 



dx 



\n(zf(a) 2 + 1) + ln(z/(6) 2 + 1) + ^ 1 _ e -2 ^ ^ + ^ 



1-e 



-2.4/fc 
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We now calculate each one of these terms individually. For the first one, we have 



k fri Jo Ja f[x) 

= 4+oi> 2 ), 



where we used Cr( — 1) — —1/12. 

For the second term, substituting si(z;x) from Eq. (l3Tj) . 



2s- 



2s- 



sm 7is 

71 

sm7rs 



7T 



fc=i 

oo 



Si(z;x 



k 



dx 



k=l 



f 2s (x)dx 



a v/l + f (x) 2 7 V/^TT 



t- s - l dt 



t(t - 4) f (x) 2 + 1 t r (x) 



l(t+l) 2 f{x) At + ll + f (x) 2 _ 



mm 

f{xf 



5s + 1 



dx + 



f(x) 2s f(x) 



a f{x)^\TfW 

b f"(x) 



1 



Qj a f(x) v /l + f'(xy 2j a (1 + /' (a;) 2 )3/ 2 



dx + 



a (1 + f(x) 2 )3/ 2 

+ 0(s 2 ). 



<ix 



For the third term, we obtain 



2s 



sin 7TS 

7T 

1 sin7rs 



2 7T 
1 



=i 

oo 

E* 



fe=l 



,-s-l 



t+T 



ln(z/(a) 2 + l) + ln(^/(6) 2 + l)" 



fi2 



[/(a) is + /(6) 



2.si 



--U2s)[f( a r+f(br} 



i 

1 

2 +S 



-2Ch(0) + 



1b/(o)+1b/(6) 



ln(2vr) 



ln/(a)+ln/(6) 



+ 0(s 2 ) 
+ 0(s 2 ), 



ZETA FUNCTION ON SURFACES OF REVOLUTION 



13 



where we used C#(0) = — |ln(27r). The fourth term involves the Euler function defined in 
Eq. 



2s- 



sm7rs 



7T 



k=l 



2s 



smrrs 



7T 



-28- 



sm7rs 



7T 

OO 



5> 2 

=1 

oo 

=1 

oo 



. k=l 



oo y - e -2 J° Si(k,z;x)dx 
Z ~'~ X 111 1 _ e -2Ak dz 

oo i _ -2 £ Si{k,z;x)dx 

z-*- 1 in — — + orr 

]_ _ e ~2Ak v ' 

z- s - 1 ln(l -e- 2Afc )^ + 0(l^ 



-2s 



. k=l 



-2 S ^ln(l-e' 2Afc ) + 0( S 2 ) 



fe=l 



= -2sln0(e- 2A ) + O(s 2 ). 
Adding these four terms gives 



C 2 (o) 



.4 



2 ' h 



-21n0(e~ 2A ) + - + ln(27r) 
6 



1 



ln/(a)+ln/(&) 



6j a f{x)^l + f'{xf 2 



(! + /' (a;) 2 ) 3 / 2 



<ix. 



ln</>(e' 



2tt 

e »• 



Applying the following formula, which will be proved as Lemma 14.3 

1 
2 

with r = A/7T, we have 

2 



— (r ) — lnr 

6 r 



C 2 (o) 



. 7T 2 , / ,n In f (a) + In f (b) 

-2 In (/>(e- 2 - + 7-r + ln(2A) + — - 1 



1 



6 £l(„\2 



f( x y 



-.dx H — 



2 

/"(*) 



6 A /(*)0^7W~~ ' 2 ^ C 1 + /'(z) 2 ) 3/2 
We now add &(0) from Eq. Q2BD with ££(0) above to obtain Eq. ()43l). 

Lemma 4.3. TTte following formula holds for the Euler function, 



dx. 



ln0(e- 2 ^)-ln0(e-^)=i 



lnr 



(47) 



(4£ 



□ 



(49) 



Proof. Consider the special case of / = R/ir on an interval of length L. Then the eigenvalues 
are 

r) + (t) ■ 



Due to the symmetry between R and L in the spectrum, <^2{s)(R, L) = £z{s)(L, R). There- 
fore, by Eq. flS}, 



C 2 (0) = -21n0 (e- 2 ^) + £| + ln(2J2) = -21n0 (e" 2 ^) + ^ + ln2L. 
Letting r = L/R, we obtain Eq. (jl9"|) . 



□ 
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5. Discussion and conclusion 



The geometric interpretation of the zeta function's residues and values can be found from 
the heat kernel expansion of the Laplacian on the surface of revolution; Refs. S3] 
gave extensive reviews of the heat kernel expansion. On a compact manifold, the heat kernel 
of the Laplacian, integrated over space, is defined as 

o(t)= Yl e ~ tXk - 

A fe eSpec(A) 

For a two-dimensional manifold with boundary, denoted by Q, the heat kernel has the fol- 
lowing expansion, 

0(f) = — - + + C + Ci Vt + dt + ■ ■ ■ . 

f Vf 2 

The coefficients are integrations of geometric invariants over Q and its boundary dQ |6J. For 
the surface of revolution, if the Dirichlet boundary condition is taken, 

M r - \ m \ r - * " 1 
4^' C ^-"87F °"6' 



Cx 



8k(a)]da, 



256^ Jm 

where \Q\ is the area of Q, \dQ\ is the total length of dfl, \ is Euler characteristic of Q, k g (a) 
is the geodesic curvature of the boundary at a, and k(a) is the Gaussian curvature of the 
surface at a. For the surface of revolution, x — 0. 
The zeta function is related to the heat kernel by 

COO = tttt / dte(t)t s -\ 

Substituting the heat kernel expansion into the equation above, we have 



ResC(l) 



ResC( 



ResC(- 



2' 

C(o) 



C_i = 




C7 i 

2 






8tt 


C = 


| = o, 



1 

47T 



2vr/(x)v / l + / / (x) 2 (ix, 
2nf(a) + 2tt/(6) 



Ci 



2v^F 
2nf(a) 

5127T 

2tt/(6) 

5127T 



5127T 



8tt 



[JfeJ(o-) - 8£;(a)]cicx 



90 



/a(o)(l + /«(o)) 



+ 



/"(«) 



/(o)(l + /«(o))a 

/"(&) 



L/ 2 (6)(l + f 2 (6)) /(6)(l + f 2 (6)) 2 . 

They agree with Theorem 14.11 completely. It confirms the validity of the contour integral 
method extended to the surface of revolution. On the other hand, the determinant of the 
Laplacian, calculated by the contour integral method, cannot be derived from the heat kernel 
expansion directly. 

Our results can also be benchmarked against an exactly solvable model with continuous 
spectrum. Consider the Laplacian on an infinitely long planar stripe of width L. Given the 
infinite volume of this configuration we need to introduce suitable zeta function densities. 
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To that end, we compactify the infinitely long planar stripe temporarily to a a finite stripe 
of width 2ttR and impose periodic boundary conditions along that direction. With Dirichlet 
boundary condition along the remaining direction, the spectrum is 

/ k \ 2 /wr\ 2 

k ' n= VR) ' G ' n = 1 ' 2 ' 3 ---- 



The associated zeta function is 

oo oo 



CM - E E 



k\ 2 /rwrN 2 



R + lT 



n=l k=— oc 

As i? — 7- oo, which is the transition to the Riemann integral, the relevant zeta function 
density is 



CM 



1 °° poo 

lim -C(s) = Y / \~°dk 

i>oo / t \ 2s— 1 

2Ch(2s-1) / (x 2 + l)-cfc - 



2s-l 



1\ / T \ 2s-l 



In particular, we have 



ResC c (l) = |, ResC Q) = Cc(0) = 0, Res^ f-^ 



and 



7T 



The stripe can be regarded as a surface of revolution with f(x) = i?, and R — > oo. For 
/(x) = i?, Theorem 14. II gives 

ResC(l) = ^, Re <Q)=-f' C(0) = 0, ResC (-0=0, 

and 

n 2 R 



C'(0) = -21n0 (V 



^ ) + 



6L 

One easily verifies, as the above computation shows, that these properties are recovered from 

CcOO = lim ^ 

As a conclusion, we have extended the contour integral method to the surface of revolu- 
tion. By the WKB expansion, we calculated the zeta function at several important points. 
The results agree with those obtained from the heat kernel expansion. We also computed 
the determinant of the Laplacian on such a surface. The WKB expansion allows the eval- 
uation of the zeta function at more points, however the expressions get increasingly more 
complicated. In the future, we can change the Dirichlet boundary condition to, say, Robin 
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boundary condition, and calculate the associated zeta function. We can also relax the con- 
dition that f(x) > by allowing f(a) = and/or f(b) = 0. This should allow the study of 
manifolds with cusps. 
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